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ABSTRACT 



Context. Most of the mass in protoplanetary disks is in the form of gas. The study of the gas and its diagnostics is of 
fundamental importance in order to achieve a detailed description of the thermal and chemical structure of the disk. 
■ The radiation from the central star (from optical to X-ray wavelengths) and viscous accretion are the main source of 

energy and dominates the disk physics and chemistry in its early stages. This is the environment in which the first 
phases of planet formation will proceed. 

Aims. We investigate how stellar and disk parameters impact the fine-structure cooling lines [Ne n] , [Ar n] , [O i] , [C n] and 
H2O rotational lines in the disk. These lines are potentially powerful diagnostics of the disk structure and their modelling 
permits a thorough interpretation of the observations carried out with instrumental facilities such as Spitzer and 
1 Herschel. 

Methods. Following Aresu et al. (2011), we computed a grid of 240 disk models, in which the X-ray luminosity, UV- 
excess luminosity, minimum dust grain size, dust size distribution power law and surface density distribution power 
law, are systematically varied. We solve self-consistently for the disk vertical hydrostatic structure in every model and 
apply detailed line radiative transfer to calculate line fluxes and profiles for a series of well known mid- and far-infrared 
cooling lines. 

Results. The [Ol]63 ^inline flux increases with increasing Lpuvwhen Lx < 10 30 erg s _1 , and with increasing X-ray 
luminosity when Lx > 10 30 erg s _1 . [Cn] 157 /xmis mainly driven by Z/Fuvvia C + production, X-rays affect the line 
flux to a lesser extent. [Nell] 12.8 /jm correlates with X-rays; the line profile emitted from the disk atmosphere shows 
a double-peaked component, caused by emission in the static disk atmosphere, next to a high velocity double-peaked 
component, caused by emission in the very inner rim. Water transitions, depending on the disk region they arise from, 
show different slopes in the correlation with the [Ol] 63 ^mline. 

0^ ■ Key words, protoplanetary disks: X-rays - disk structure - IR fine structure line emission 

in 

p; 

OV 1. Introduction component. They found no accretion rate higher than 10~ n 

M yr _1 for objects older than 10 Myr; assuming an expo- 

CN ■ Protoplanetary disks are the intermediate step between the nential decay they infer a timescale of 2.3 Myr. Mechanisms 

^ I cloud collapse and the planetary system stage. The under- such as photo-evaporation can r emove the gaseous disk 

J> . standing of the chemical and physical properties of disks on short tim e scales (~ 10 5 yr~*) (iHollenbach et all Il993t 

'O ! P rovides tne initial conditions for planet formation. The lHaisch et ail 1200 li lAlexandeii 120081 ). The accretion rates 

r> ! study of the gaseous component has gained increasing in- found bv lFedele et all (|2010l) would then indicate that the 

^ 1 terest in the last few years as it dominates the mass budget gas and dust lifetime are very similar. 

in the disk in its early stages and because the increased The gas densities and temperatures in disks are suitable 

sensitivity of new instrumentation allowed us to observe it, f or the excitation of IR fine-structure transitions of species 

e.g. Spitzer, Herschel, VLT, ALMA. suc h as O, C+, Ne+ and Ar+. Different molecular species 

Observations with the Spitzer Space Observatory were (e.g. CO, H 2 0) probe different disk regions, thereby allow- 

used to estimate that the lifetime of the inner "dusty" disk ing to build a coherent picture o f the whole d isk, from a 

is of t he order of 10 Myr (IStrom et all . Il989t lHaisch et all , complete suite of observation fsee lBerginl (|2009l ) for a com- 

12001k iHernandez et all 120081) . Finding the corresponding prehensive review). These are among the diagnostic tools 

timescale for the gaseous component of the disk is rather that are used to infer disk and stellar properties from ob- 

complex, because there is not a single gas tracer that can be servations with IR satellites like Spitzer and Herschel. The 

used; instead recent modelling efforts suggest that a suite of latter covers some major disk cooling lines such as [0 1] 63 

several gas tracers may be necessary to cover th e wide range A*ni and [C 11] 157 /im. 

of chemical and excitation conditions in disks (jGorti et all The origin of [0 1] 63 /mi emission from protoplane- 

l2011tlKamp et all 1201 1[ ). [Fedele et al.l (|2010D looked at ac- tary disks is not unique: it can originate from excita- 

cretion signatures to indirectly study the lifetime of the gas tion in the disk and from the interaction of jet /outflows 
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with the circumste llar environment (Ccccarc lli et allll997t 
iNisini et all Il999t) . Without spatially resolved observa- 
tions, the contr ibution of these tw o mechanisms is difficult 
to disentangle (jPodio et all l2012h . 

[O i] 63 /im is recognized to be a temperature tracer for 
the surfa ce layers of protopla netary disks between 10 and 
100 AU (IWoitke et all 12009ft. An exp loratory study with 
the DENT grid ()Woitke et all (|2010l) . 300000 disk mod- 
els based on the junction of the thermo-chemical code 
ProDiMo and the radiativ e tran sfer code MCFOST) was 
carried out by iPinte et all (|2010f ) for a sample of T Tauri 
and Herbig stars. They show that FUV radiation, poten- 
tially arising from accr etion onto the st a r, play s a dominant 
role in T Tauri stars. Me iierink et al.l (|2008f) investigated 
the effect of X-rays, finding a correlation between the line 
flux and the X-ray luminosity. iGorti fc Hollenbachl 



calculated the [O i] line for models with different FUV and 
X-ray luminosity, finding a major impact on the line flux by 
I/puv- lAresu et al.l (|201ll ) considered X-rays and FUV ra- 
diation, finding that for their fiducial solar-type star model 
X- rays become imp or tant a bove Lx ~ 10 30 erg s _1 . 

iGlassgold et all (l2007h , f ollowed by work from 
IGorti fc Hollenbachl (|2008l) and lErcolano fc Owed (|2010D . 
predicted emission of ionized neon at 12.8 /im in the 
atmosphere of an X-ray irradiated protoplanetary disk as 
a probe of physical conditions of the hot atmosphere (~ 
4000 K) in the inner disk (< 20 AU). The line was detected 
for the first time the same year wi th the infra-red spec- 
trometer (IRS) on b oard of Spitzer dPascucci et all l2007t 
lEspaillat et all 120071; iLahuis et all 12007)1 The origin of the 
[Neil] 12.8 /zm line is still under debate. iGiidel et all (|2010l ) 
analysed a collection of more than 50 [Ne n] detections, and 
suggested different mechanisms responsible: emission from 
(a) the disk, (b) the photo-e vaporative flow from the disk 
surface and (c) the j e t itsel f (IShu et all 119941: [Shane et al.l 



; and (c) tne l et itseit 
ISacco et all (|2012h 



and Baldovin-Saavedra et alj 



analysed high spectral resolution [Neil] lines. They 
infer from the position of the peak with respect to the 
stellar velocity and from the line profiles that the lines are 
mostly excited both in photoevap orative flows and in jets. 
Spectrally resolved observations (jNaiita et all I2009D also 
found a few sources where the [Ne n] is most likely emitted 
in a disk atmosphere. 

The study of the gaseous disk remains an important 
challenge, where modelling can provide further insights. 
Thorough computing of the disk thermal and chemical 
structure, when coupled to observations of gas tracers of 
different disk regions, can be used as a powerful baseline 
to predict the fate of protoplanetary disks. Many proto- 
planetary disk models include X-rays, as there is general 
consensus that this radiation plays an important role in 
the chemistry, and in the thermal balance of the disk. 

In this context several groups developed numerical 
codes dedicated to the mode ll ing of protoplanetary 
disk (IGorti fc Hollenbachl 12004 IGlassgold et all [2 007: 



Nomura et all. l2007t lErcolano et all. 120081: IWoitke et al 
2009 : IGorti et all 120091: iHollenbach fc Gortil. 12009: 
Woods fc Willacvl. 12009ft . IWoitke et all (|2010h and 
Kamp et al.l < 201lh investigated the effect of a wide 



parameter space on disk properties and studied the line 
emissi on of species such as CO, O and C + . iMeiierink et al.l 
(|2008tl studied the effect of X-rays on Ne+ and Ne 2+ , 
different ionization stages of S, C, C + as we ll as O for a 
low mass star disk model. IGorti et al.l (|2009h investigated 



the role of FUV (6 < hu < 13.6 eV), EUV (13.6 eV < 
hv < 100 eV) and X-rays (hv > 100 eV) on the emission 
of [Aril] at 7 ^m, [Nell] at 12.8 /im, [SI] at 25 /im, [Fell] 
at 26 /im, [OI] 63 fim and also pure rotational lines of 
H 2 and CO. They also investigated the properties of 
EUV/X-ray driven photo-evaporative flows in the upper 
layers of protoplanetary disks and their potential impact 
on diagnostic lines. 

iFogel et all (|201lD studied the effects of grain settling 
and Lya line scattering on the chemistry finding that both 
effects arc important for molecules like CO, CN, HCN and 
H2O, and in general for the carbon and oxygen molecu- 
lar chemistry. iHeinzeller et al.l (|201ll ) showed that viscous 
heating and turbulent mixing also impact the molecular 
layer, thereby modifying the thermo-chemical conditions. 

The scope of t his paper is to follow up earlier work of 
lAresu et al.l ([201 ll ) and to study the relative importance of 
X-rays and FUV in setting the thermal and chemical con- 
ditions in the disk. This work will investigate how selected 
atomic and molecular diagnostics respond to the differ- 
ent energy input. Meijerink et al. (2012), hereafter Paper I, 
show the respective roles of X-rays and FUV stellar radia- 
tion on the disk physical and chemical structure. The grid 
consists of 240 models, where we varied Lx, £fuv, dust 
properties (a m i n , a pow ) and the power law for the surface 
density distribution (e). We provide, here in the second pa- 
per, an extensive study of the line diagnostics such as [0 1] , 
[C 11] , [Ne 11] , [Ar 11] and water and investigate the thermo- 
chemical conditions behind the excitation of these lines. 

This paper is structured as follows: in Sect. 2 we briefly 
recall the setup of our model grid. In Sect. 3 we list and 
explain the results obtained, in Sect. 4 we compare our 
findings to those of previous works. In Sect. 5 we summarize 
and provide an outlook for future work, in Sect. 6 we draw 
the conclusions. 



2. The grid of disk models 

We use the thermo-chemical code ProDiMo to compute 
thermo-chemical disk models in hydrostatic equilibrium. 
The parameters adopted are listed in Table [TJ We consider 
a Sun-like T Tauri star surrounded by a 10 -2 M Q disk. To 
study the impact of the stellar radiation on the disk physics 
and chemistry we vary X-ray (Lx) and FUV (Lfuv) lumi- 
nosities, dust properties (minimum dust size and dust size 
distribution power law) and the surface density distribu- 
tion power law. The dust to gas mass ratio is kept fixed 
(p=p d /p g =0.01). 

Every disk model is composed of 100 radial and vertical 
points where we solve for the thermal and chemical balance. 
The resulting vertical temperature profile is then used to 
compute the vertical disk structure until hydrostatic equi- 
librium is reached (see Woitke et al. 2009). A self consistent 
treatment of the disk flaring is needed to compute reliable 
optical depths through the disk. 

The source of the FUV and X-ray radiation is assumed 
to be centred at the stellar position (a representative spec- 
trum is shown in paper I). No X-ray scattering was con- 
sidered. We ran the grid on the Millipede cluster of the 
University of Groningen, see paper I for more information. 
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Tabic 1: Parameters used in the models. 



Quantity 


Symbol 


Value 




M* 


1 Mr-, 


ILllcCLl Vc LcIIljJcI cLL Lilt; 


-t cff 


t 1 U IV 


O Lclldl ILlXIllXiUolLy 


T 


1 T r~. 


J * loJ\ AXJLCU3D 


Ma- 1 


o ni Mr^ 

U.Ul 1 V± [V) 


X-ray luminosity (0.1-50 keV) 


T „ 


U, 1U , 1U 






iu , iu erg/ s 


r V V luminosity 


r 

-Lfuv 


-1 n 29 1 n 30 

10 ,10 , 






10 31 , 10 32 erg/s 


Inner disk radius 


Tin 


0.5 AU 


Outer disk radius 


r ut 


trnn ATT 

OOO AU 


Surface density power law index 


£ 


l.Uj 1.0 


Dust-to-gas mass ratio 


pd/p 


U.Ul 


Nlin. dust particle size 


f^min 


U.l, U.O, l.U fJJJX 


Nlax. dust particle size 


Imax 


IU UlU 


Dust size distribution 


^pow 


2.5, 3.5 


power index 






Dust material mass density 


Pgr 


2.5 g cm -3 


Strength of incident ISM FUV 


X ISM 


1 


Cosmic ray ionization rate of H2 


CCR 


5 x 10" 17 s _1 


Abundance of PAHs relative 


/pah 


0.12 


to ISM 






viscosity parameter 


a 






2.1. Line radiative transfer 

The line rad i ative transfer is calculated as explained in 
IWoitke et al.l (|2011h . section A. 7. We chose an inclination 
of 45 degrees and a maximum range of ±30 km/s for the 
width of the lines. The line profile is calculated taking into 
account thermal and turbulent broadening, the latter is set 
to uturb =0.15 km/s. The flux is calculated for a disk at 140 
pc with an inclination (i) of 45°. For a subset of 20 models 
with varying Lxand Lfuv, and a m i n = 0.1, a pow — 3.5 and 
e= 1.5, the line fluxes have been calculated also for i = 0°, 
60°, 75°, 90°. The transitions described in this paper are 
listed in Tabled 



2.2. Data 

The chemical network counts ~1500 reactions, part of it is 
composed of UMIST reactions ()Woodall et all 120071) . while 
the reactions involving X-ray primary and secondary ioniza- 
tion of all the elements are taken from lMeiierink fc Spaansl 
(|2005l) and lAdamkovics et al.l f|2011h . The treatment of 
molecular di ssociation due to X -ray absorption is explained 
in Table 1 in lAresu et al.l (|2011[ ). Charge exchange reactions 
between single or double ionized ions with neutral species 
are taken from lAdamkovics et al.l (|201l[ ). 

In Table [5] we show the references for the collisional 
rates used for oxygen, ionized carbon, ionized sulfur and 
water. The electron excitation cross secti on for Ar + and 
Ar 2+ were computed in the Iron Project (jHummer et all 
Il993l) . Collisional rates for Ne + and Ne 2+ are taken from 
the CHIANTI database [Dere et al.l (fl997l . [2009]) . 



3. Results 

In this section we describe the resulting line fluxes for 
[Ol]63 /xm, [C11] 157 pun, [Neil] 12.8 /imand [Aril] 7 pun. 
For a better understanding of the thermal and chemical 



conditions under which these species emit, we also give a 
description of the mass averaged gas temperature (defined 
below) and the total species mass. In addition, ProDiMo 
computes both the location of the line emitting region (de- 
fined below) and the FWHM of the line. These quantities 
will be used to investigate how the different stellar and disk 
parameters affect the emission of the coolants. 

The dominant factors controlling these lines are £fuv 
and Lx- Hence the quantities we will describe in the fol- 
lowing are averaged among models with fixed Lpuv and 
L x - The effect of varying a m i n , flpow and eis then shown in 
form of an error bar. Defining a "series" as a suite of mod- 
els with fixed Lxand LFUvbut different a m - m , a pow and e, 
we count 12 models in each series. This way, we consider 
for a given physical quantity Q in the series k, the simple 
arithmetic mean Q k . To avoid overshoot in the error bars 
due to a handful of models (n < 10) that did not reach opti- 
mal global convergence, we make use of the weighted mean 
for plotting purposes. The squared distance of qk,i from the 
arithmetic mean is then calculated and the inverse is used 
as weight (uJk.i — , ; ) to compute the weighted average 
of the quantity Q for the k-th series: 



l^i=i Ik,. 



E12 



We computed the variance a\ for each series and the error 
bars in the plots represent then the 2<r deviation from the 
Qw- This approach is used for all the quantities discussed 
in this paper (flux, mass averaged temperature, mass of 
species sp, R m etc.) 

We also use the mass averaged temperature for a species 
sp, defined as: 



(T) 



sp 



J T gas (r, z) n sp (r, z)m sp dV 
J n sp (r, z)m sp dV 



(2) 



We calculate the mass averaged temperature only in the 
emitting region of the considered species: applying a ID es- 
cape probability line radiative transfer, we calculate the 
radial (i?; n and i? ou t) and vertical coordinates that enclose 
the region where up to half of the total flux of the line is 
emitted. This volume is then used to perform the integra- 
tion in Eq. |3J This allows to relate the temperature directly 
to the flux of the emitted line. 

3.1. Oxygen fine-structure emission at 63 pm 

The lower left panel of Fig. [T] shows how the line flux 
changes as a function of the X-ray luminosity; color coded 
is the UV luminosity. Plotted is the weighted flux for mod- 
els with a given X-ray and FUV luminosity and the error 
bars represent the 2cr scatter due to the variation of the 
flux with respect of the other parameters in the grid. 

In the absence of X-rays, [0 1] emission is driven by 
Lfuv'- the 63 pim line flux increases by a factor ^10 
from the lowest FUV model with Lpuv = 10 29 erg s _1 to 
Lfxjv = 10 31 erg s _1 and another factor ~ 10 when Lfuv is 
increased to 10 32 erg s _1 . For constant Lfuv, X-rays always 
start to impact the [0 1] 63 pim emission beyond a thresh- 
old of Lx~ 10 30 erg s . The only exception is the highest 
LfuVj where [Oi]63 pm stays constant up to Lx^ 10 31 . 
For the highest X-ray luminosity, [0 1] reaches a plateau, 
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Table 2: List of transitions considered in the line radiative transfer. From left to right we list the species name, the 
transition label, the excitation temperature, the Einstein coefficient for spontaneous emission, the collisional partners 
included in our model and approximately the region of the disk that the transition traces. 



Species 


Transition 


A [/xm] 


T CX [K] 


1 7 T\ 

A (s- 1 ) 


Coll. partner 


Region 


Ref. 


O 


J Pi-^P 2 


63 


230 


8.91xl0 _b 


H a , H, H + , e" 


20-200 AU 


1,2,3,4 


O 


J Po- J Pi 


145 


326 


1.75xl0~ b 


Ha, H, e" 


20-200 AU 


1,2,3,4 


c+ 


ri/2" *3/2 


157 


91 


2.30xl0~ b 


Ha, H, e" 


Outer Disk 


1,2 


Ne+ 


2-p 2-p 

ri/2- ^3/2 


12.8 


1122 


8.59XRT 3 


H, e" 


Inner Disk 


5,6 


Ne y + 


3 Pl- 3 P 2 


15.5 


925 


5.97xl0~ a 


H 


Inner Disk 


5,6 


Ar+ 


2-p 2-p 

ri/2" -T3/2 


6.9 


2060 


5.30x10^ 


H 


Inner Disk 


7 


o-H 2 


2l2 — > loi 


179 


114 


5.59x10^ 


H 2 , H 


Outer disk 


1,8 


0-H2O 


818 — > 7o7 


63.3 


1071 


1.751 


H 2 , H 


Hot belt 


1,8 


0-H2O 


845 — > 7l6 


23.9 


1615 


1.05 


Ha, H 


Inner wall 


1,8 



References for the collisional rates: 



,2()0(ili . (4) 

(jStorzer fc Hollenbachl . l2000t ). (5,6) iDere et all (|1997l |2009| ).~(T) iHummer et alT(|1993l) . (8) The rates for collisions with 



atomic hydrogen are a scaled version of the H2 data. 



-18 ' — 1 



1 



-uv 
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-uv 



1e29 erg/s 
1e30 erg/s 
1e31 erg/s 
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3 
O 
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25 
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L uv = ie29 erg/s 
L uv = 1e30 erg/s 
L uv = 1e31 erg/s 
1e32 erg/s 



i 

_i 



uv 



29 



30 

L x [erg s" 1 ] 



31 



32 



Fig. 1: From lower to upper panel on the left: flux of the [Ol] 63 /im line, oxygen mass averaged temperature and total 
oxygen mass in the disk versus X-ray luminosity. See Sect. 3 for the description of the error bars. Since the line is 
optically thick, its flux is sensitive to the thermal conditions in the disk. From lower to upper panel on the right: [Ol] 63 
/im emitting region inner and outer radius and FWHM of the line. As the emitting region gets pushed to larger radii, 
the FWHM of the line decreases accordingly. 
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Fig. 2: Left panel: O abundance (relative to hydrogen nuclei). The black contours are iso-temperature curves at 200 and 
1000 K. Central panel: optical depth of the line (blue) and of the continuum (black) in the top mini-panel, below is the 
cumulative flux of the line, which shows which percent of the final flux is built and where, the lower mini-panel shows the 
oxygen density, the contour indicates where the contribution to the line flux goes vertically from 15% to 85% at every 
radius. Right panel: FUV (top) and X-ray (bottom) heating rates. These rates are obtained integrating the volumetric 
heating rate T over the volume of the [OI] emitting region for every model. Averages and error bars have the same 
meaning as described in section [3] The first two panels are taken from a representative model with Lx= 10 30 erg s , 



10 erg s 



0.1 fim, a pow = 3.5 and e= 1.5. 



which represents the maximum value of ~ 10~ 15 W/m 2 in 
our series of models. 

The middle left panel of Fig. [JJ shows the mass averaged 
temperature of atomic oxygen vs Lx in the region of the 
disk where 50% of the total flux is emitted. The behaviour 
is similar to the one discussed for the line flux: the tem- 
perature increases with increasing input energy (see also 
Fig. 2 in paper I). In the UV only models (T) Q slightly 
increases from T = 35 K in the low Lpuv models, until it 
reaches about 80 K for Lfuv= 10 32 erg s _1 . The models 
with the highest FUV luminosity show also the bigger scat- 
ter in temperature (AT± 20 K). When we let Lx increase, 
the temperature for the low Lpuv models starts to increase 
when Lx^ 10 30 erg s -1 . In higher Lx models, the tem- 
perature in the [O i] emitting region is set by X-ray energy 
deposition, regardless of the Lpuv contribution. 

In the upper left panel of Fig. [TJ we show the variation 
of the total oxygen mass in the disk. The OI mass only 
varies within a factor 2.5 across all the models. The scatter 
caused by different dust properties (a m i n , a pow ) and mass 
distribution in the disk (e) for a series of 12 models is at 
most ~ 2.3, hence is higher than the variation of the av- 
eraged oxygen mass with Lpuv m different series of models 
(maximum variation of a factor 1.4). On the other hand the 
X-ray luminosity does not affect the total oxygen mass in 
the disk. This is also explained in section 4.5 of paper I. 

The [O i] emission is optically thick in all models. Its 
flux is then sensitive to temperature variations. The mutual 
FUV/X-ray contribution in creating the thermal conditions 
that drive [O i] emission becomes clear from Fig. [2] (right 
panel) . We calculated the volume integral for the volumetric 
heating rate of X-rays and FUV, obtaining the respective 
heating rates in the oxygen emitting region: 

*i = y" T i {r,z)dV [ergs" 1 ] (3) 

Here Tj is the volumetric heating rate in the i-th model for 
X-rays or FUV radiation in units of erg s _1 cm -3 . We used 



our previous definition of the oxygen emitting region as 
integration volume. The FUV emitting rate is the sum of 
all relevant FUV related processes: photoelectric heating, 
PAH heating, carbon ionization heating and H2 dissocia- 
tion heating (see Woitke et al. 2009). In the lower panel, 
we show how the X-ray heating rate scales roughly pro- 
portional to Lx- ^x-rays ~ 0.01 Lx- In the upper panel, 
we show the FUV heating rate, we find approximately that 
^fuv ~ 0.01 — 0.001 Lpuv- Note how for low Lpuv models 
increasing X-rays causes progressive decreasing in the FUV 
heating rate. This is due to a general decrease of the local 
FUV radiation field at high Lx, which causes FUV heat- 
ing to decrease, especially PAH heating, which is the main 
FUV related heatin g process in the ox ygen emitting re- 
gion (see appendix in lWoitke et al.l ((201 lh for the treatment 
of PAH heating in ProDiMo). X-ray energy deposition be- 
comes comparable with FUV when Lx ~ 10 30 erg s _1 . It is 
the competition between PAH and Coulomb heating that 
sets the temperature in the oxygen emitting region, thereby 
impacting directly on the flux of the line. 

The location of the [O i] emitting region varies through 
the models, especially for the highest values of Lx- Inner 
and outer radius of the emitting region generally both in- 
crease with increasing X-ray luminosity (lower and middle 
right panel of Fig. [JJ). Inner and outer radius also increase 
with Lfuv- Thereby suggesting a dependence on the total 
energy input (in the X-rays and FUV band). 

In the top right panel of Fig. [TJ we show the full width 
half maximum (FWHM) of [O i] throughout the model pa- 
rameter space. The FWHM decreases for increasing Lx, 
which is in accordance with the changes in the line forming 
region. 

3.2. Carbon 

The ionized carbon fine-structure line flux at 157 /im 
clearly depends on Lpuv for the UV only models. UV ra- 
diation (A < 1100 A) sets C + abundances in the outer 
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Fig. 3: Lower and central panel: inner and outer radius of the emitting region. Top panel: full width half maximum of 
the line. See Sect. 3 for the description of color coding and error bars. LEFT figure: [O i] at 63 /xm. The oxygen emitting 
region is generally pushed further as Lx increases, this is also reflected by the FWHM, which on average slightly decreases 
with Lx- RIGHT panel: [Cn]at 157.7 /xm. Ionized carbon traces the outer disk, its emitting region extends inward as 
Lpuv increases. X-rays do not impact considerably Rm, i?ouT or the FHWM of the line. 



disk (see paper I). The increasing X-rays can contribute to 
the line intensity for the lowest two FUV luminosities for 
Lx < 10 30 erg s _1 . After this value a "plateau" is reached: 
even though X-rays increase the temperature in the ionized 
carbon emitting region, the column density in the same re- 
gion of the disk is decreasing (see the first two panels in the 
last three rows of Fig. A. 24 in paper I), hence the latter ef- 
fect compensates for the temperature increase. For the high 
FUV models the column densities do not vary for models 
with different Lx- Hence the line flux is mainly increas- 
ing with LpuVj with a minor contribution of Lx, through 
thermal effects, as in the [Oi] 63 /imcase. 

For Lx< 10 30 erg s^ 1 , the C + temperature is UV driven 
and increases for the highest -Lfuv(10 30 erg s _1 ) by a factor 
2 compared to the low Lpuv (10 29 erg s -1 ). At higher Lx, 
the C + temperature is entirely controlled by X-rays driving 
it up to ~ 100 K in the most extreme case. 

In the top left panel of Fig. [31 we plot the total mass of 
C + as a function of X-ray luminosity. The role of Lfuv in 



the formation of ionized carbon is clearly shown: the total 
mass spans roughly one order of magnitude from 3.2 xlO -8 
M Q for the lowest Lfuv to 3.2 xl0~ 7 M Q for the highest 
one. X-rays and the other free parameters affect the total 
C + mass to a lesser extent: for a given LpuVi mc+ changes 
at most by a factor 2. 

The [C ii] emission is generally very close to being opti- 
cally thick. In this regime the line is sensitive both to the 
column density of the species and to the gas temperature 
(see paper I, A. 24). Indeed in our models, the ionized car- 
bon flux is mainly dependent on its total mass in the disk. 
However, if the temperature of the emitting region changes, 
the flux varies accordingly. These two effects arc caused 
from a combination of FUV radiation, which controls the 
production of C + (top left panel of Fig. [3j) , and from X-ray 
radiation, which, similarly to oxygen, contributes to the 
thermal balance in the ionized carbon emitting region. 

The location of the emitting region is mainly dependent 
on the FUV flux (Fig. [3]). In the high luminosity models, C + 
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Fig. 4: From lower to upper panel: flux of the line, mass averaged temperature and total species mass in the disk versus 
X-ray luminosity. See section 3 for the description of the error bars. LEFT figure: ionized neon fine-structure line emission 
at 12.8 /Ltm. The line is produced in a high temperature X-ray heated environment, where the density is low (nn ~ 10 6-7 
cm -3 ). The line is optically thin and hence very sensitive to the total ionized neon mass, which correlates with Lx- 
Hence, in our models, [Ne n] correlates with Lx- RIGHT figure: ionized argon fine-structure line emission at 6.9 /jm. 
[Ar ii] behave very similarly to [Ne n] , and it can be used as a tracer for the same region. 



emits closer to the central star, starting from r ~ 80 AU. 
For constant LfuV; the inner radius of the emitting region 
i?i n , increases only weakly with Lx- However, it shows a 
strong inverse dependence on Lpuv : it increases by a factor 
3 from high to low Lpuv models. The outer radius of the 
emitting region shows the same trend. The different radial 
location of the ionized carbon emitting region causes the 
FWHM of the line to change as shown in the upper right 
panel of Fig. [3] The variation is less than a factor 2 from 
high to low Lpuy models. As expected, the FWHM is not 
noticeably affected by Lx- 

3.3. Neon 

Ionized neon is expected to be sensitive to the X-ray ra- 
diation, because neutral neon has a primary ionization for 
the K-shell of ~ 900 eV. In the lower left panel of Fig. gj 
we show the flux of [Neil] versus Lx- UV only models are 
not shown as they do not produce significant [Neil] emis- 
sion. The dependency of [Ne n] on parameters other than 
Lxis weak, as reflected in the small error bars for every 
series of models. The linear fit gives a slope of 0.81 for the 
correlation of the line flux with Lx- 

In the upper panel of Fig. HJ we show the total Ne + 
mass in the disk. The higher Lx, the more ionized neon is 
produced. Other parameters also affect the total Ne + mass 
budget, but the overall trend is driven by X-rays. The linear 
fit gives a slope of 0.75. 

The lower and upper left panel of Fig. [1] shows clearly 
that ionized neon fine-structure line emission is controlled 
by X-ray ionization of neon. The mass averaged temper- 



ature, not shown, in the [Ne n] emitting region ranges be- 
tween 2000 K and 5000 K. The electron fractions in the 
upper layers within 20-30 AU are generally high (x c \ < 0.1) 
but the total electron density where most of the line forms 
is below n cr (few times 10 5 cm -3 ) causing the line emission 
to be mostly not in LTE and optically thin. Its intensity 
is then regulated by the column density of the species at a 
given radius, which is clearly set by the X-ray luminosity 
(see paper I, A. 27). 

The lower and central left panel of Fig. [5] show how 
[Ne n] emission is confined within 20 AU. The emitting re- 
gion is pushed further out as Lx increases. As a result, the 
FWHM of the line decreases accordingly (upper left panel, 
fig. from 25 to 12 km/s. 

Collisions with H 

Collisions with H are potentially important in the exci- 
tation of the fine-structure lines of Ne + . In our grid, we 
only considered collisional excitation of ionized neon with 
electrons (Ccoii.H = Ci = 0). A c ollisional excitat i on rat e 
for atomic hydrogen is given in iBahcall fc Woll (Il968h : 
(coii.H ~2x 10~ 9 cm 3 s _1 (£2)- This value does not take 
into account the spin interaction between H and the tar- 
get electron and hen ce overestimates the r ate by almost 
a factor 10. Following iMciicrin k et al.l ()2008f ). we therefore 
adopt (coil, h ~2x 10" 10 s" 1 (£3). Using C2 and C3 we run 
a subset of three models with increasing Lx from 10 29 to 
10 31 erg s , all other parameters remain unchanged from 
what refer to as the " standard model" (paper I) . Table |3] 
shows how the collisions with H tend to increase the line 
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Fig. 5: Lower and central panel: inner and outer radius of the Ne + emitting region. Top panel: full width half maximum 
of the line. See section 3 for the description of color coding and error bars. 
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Fig. 6: Left panel: Ne + abundance (relative to hydrogen nuclei). The black contour encloses the region in the disk where 
the electron density exceeds the critical density for the collisional excitation of the 12.8 /imline. Central panel: optical 
depth of the line (blue) and of the continuum (black) in the top mini-panel, below is the cumulative flux of the line, 
which shows which percentage of the final flux is built and where, the lower mini-panel shows again the ionized neon 
density, the contour indicates where the contribution to the line flux goes vertically from 15% to 85% at every radius. 
Right panel: line profile and image for a distance of 140 pc and an inclination of 45°. 
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Table 3: Line fluxes for the 12.8 /tm transition of Ne + con- 
sidering only collisions with electrons and two different 
rates for excitation collisions with H. 



L x [erg s _i ] 
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10 31 


C [cm 3 s^ 1 ] 


Flux [W m - *] 


Flux [W m - *] 


Flux [W m - *] 


2 x io _a 


4.5e-19 


2.6e-18 


1.6e-17 


2 x icr iU 


2.3e-19 


1.3e-18 


7.4e-18 





1.6e-19 


1.0e-18 


5.7e-18 



flux with respect to the standard models, by a factor ~ 3 
for CcoiLH = C2 and by a factor 1.3-1.4 for Ccoii.H = (3- 

Collisions with hydrogen affect the total [Ne 11] flux by at 
most a factor of a few, and the trend is systematic through 
all models. Given the uncertainties in the collisional rates, 
we did not further consider excitation by H collisions in the 
grid of disk models presented here. 

3.4. Argon 

In Fig. 2] (right figure, lower panel), we show the fluxes for 
ionized argon fine-structure emission at 6.9 /mi. In the same 
figure, central and upper panel show the mass averaged 
temperature and the total Ar + mass in the disk. Except 
for the absolute value of the line flux (F^ r + ~ 10 x L Nc +), 
all the results described above for Ne + apply. The ther- 
mal properties in the [Ar 11] emitting region are the same 
as described above for Ne + . The total mass of Ar + follows 
the same behaviour, in terms of magnitude and dependence 
from Lx, as found for Ne + . 

3.5. Water 

We included in the line radiative transfer also a few rota- 
tional water transitions (Table [5J . Given the high number 
of transitions and their broad range of excitation tempera- 
tures, the water molecule emits over a wide range of radii, 
tracing disk regions that can differ vastly in terms of chem- 
ical and temperature properties. Understanding the water 
excitation conditions offers the chance to ext ract informa- 
tion o n the radial disk structure. Following IWoitke et al.l 
(2009J), we choose to divide the disk into three representa- 
tive zones: the cold water in the outer disk, the warm to hot 
water in a zone above the midplane ranging from i?j n to well 
beyond the snow line and the warm to hot main gas phase 

In Table 



water reservoir behind the inner rim (Fig. IA.1I) 
we list the transitions which are well suited to trace these re- 
gions: 179 /im (cold region), 63.3 /im (hot surface water) and 
23.9 /im (inner wall water). In Fig. [3 we show how the dif- 
ferent transitions respond to variations in Lx and Lpuv hi 
our grid. We can identify three different behaviours. The 
line emitted in the inner wall (23 /iinline, lower panel) is 
not altered by the X-ray radiation, unless Lx= 10 32 erg 
s _1 , but only responds to Lpuv- The line formed in the hot 
water layer shows an oxygen-like behaviour: X-rays increase 
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Fig. 7: From lower to upper panel: water line fluxes for the 
23.9 /im, 63.3 /imand 179 /tmtransition versus the X-ray 
luminosity. Colour coding is the same as shown in Fig. [TJ 
The X-ray impact is more evident as we move toward the 
outer disk. 



water density distribution. We overplot the emitting regions 
for the 0-H2O transitions at 179, 63 and 23 /im. The boxes 
enclose the disk region where the contribution to the to- 
tal line flux for a given transition is about 50%. All the 
lines discussed here are highly optically thick, hence they 
are sensitive to temperature variations. We also show the 
[0 1] 63 /tm transition emitting region previously described 
for comparison. Red contours are iso-temperature curves at 
Tgas = 2000, 200 and 50 K, black contours indicate the r(l 
keV) = 1 and AV = 1 lines. 



the line fluxes beyond a given threshold, 10 3 erg s (vs o -H 2 23 9 \xm 



10 30 erg s _1 for oxygen). The outer disk line flux (179 /im) 
also increases with increasing Lx, and it shows no depen- 
dency on Lpuv- To understand these trends, we show in 
Fig. I A. II the emitting regions for this lines in a subsamplc 
of models: starting from the top left panel, Lpuv increases 

erg s _1 while Lx increases 
1 . The grey scales shows the 



horizontally from 10 to 10 



V'52 



vertically from to 10 erg s 



This line is emitted in the inner disk in a high density region 
(rtH ~ 10 10 cm -3 ), where the disk faces directly the stellar 
radiation and the water abundance relative to the hydrogen 
atomic nuclei is 10 ~ 4 . The emitting region is optically thick 
to the X-ray radiation but coincides with the AV=1 line. 
Increasing Lpuv indeed seems to affect the total emitted 
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flux, while I/x radiation does not play a role (Fig. [7J lower 
panel). The temperature in the line emitting region is set 
by FUV heating. 



that of the 23 /mi line (0.3), even though the cause is 
fundamentally different. 



0-H2O 63.3 /mi 

The line traces the high temperature regions above the mid- 
plane in the inner disk. In this hot water layer, we can 
recognize a heating pattern similar to that discussed for 
oxygen in Sect. 3.f (Fig. [71 central panel). Only after a 
given luminosity threshold, 10 31 erg s _1 , X-rays are able to 
contribute to the heating in this region; below Lx= 10 31 
erg s" 

region that lays between r(l keV) = 1 and the AV = 1. For 
X-rays to contribute to the thermal balance in the hot water 
belt, Lx has to be higher than 10 31 erg s" 1 and Lx >Lfuv- 



4. Discussion 

In this section, we discuss our results, comparing them with 
previous works and analysing their diagnostic potential in 
the interpretation of disk observations. 



4.1. [Oi] at 63 /mi 



FUV heating dominates. The line is formed in the The [° ^ flux is both dependent on L x and L F uv- Except 



o-H 2 179 urn 

This transition is excited in the outer disk (r ~ 30-200 
AU). The line flux is more sensitive to the X-ray radiation 
because water forms in the outer disk through ion-molecule 
chemistry (see paper I) and X-rays dominate the thermal 
balance there. The excitation conditions for this transition 
are then shifted to higher layers with respect to the FUV 
only models; these are layers where the X-ray optical depth 
approaches one at 1 keV (and AV < 1) and the temperature 
is set by Coulomb heating. Fig. lA.II shows how the emitting 
region "moves" upwards as Lx increases. Note that apart 
from the extreme Lx model, the 179 /mi flux only increases 
by a factor ~ 3. 



Correlation with [Oi] 63 /mi. 

In Fig. [5J we show the water lines fluxes versus the [0 1] 63 
/mi line flux. At low X-ray luminosities, the 23 /online 
weakly correlates with [0 1] 63 /mi. The correlation is due 
to the increase of Lfuv, which heats both emitting re- 
gions, causing both line fluxes to increase. The slope of the 
correlation would be lower than one because, for a given 
change in LpuVj the [Oi] flux increases more than the wa- 
ter 23 /mi flux, since the latter is shielded by the inner rim. 
At higher X-ray luminosity, the line flux for [0 1] increases, 
while the 0-H2O line does not change anymore. Increasing 
the FUV luminosity no longer affects the oxygen flux, while 
the 23 /tm water line still increases with Lfuv- The X-ray 
models then fall below the previous correlation, causing 
more scatter in the plot (Fig. [51 left panel). Although we do 
not provide a " best fit" , we still show a line for the whole set 
of models, to illustrate how the slope would change when 
X-ray models are included. 

The same effect can be seen in the correlation between 
water emission at 63 /mi and [0 1] 63 /mi. Except in this 
case, the correlation at low X-ray luminosity is almost lin- 
ear (0.96), because an increasing FUV luminosity causes a 
uniform increase in both line fluxes. For Lx~ 10 30 erg s _1 , 
the global correlation is again shallower (0.74). 

No correlation is present between 0-H2O 179 
/mi emission and [Oi]at low Lx- The former line is 
indeed not affected by Lfuv- Though, when Lxis higher 
than ~ 10 30 erg s _1 , both line fluxes increase. The slope 
of the correlation is small and ends up comparable with 



in the presence of strong FUV excess (Lfuv= 10 er g s 
our results suggest that Lx has to be larger than 10 30 erg 
s _1 for the [Oi]flux to increase notably. At the threshold, 
Coulomb heating contributes as much as FUV heating to 
the thermal balance in the line emitting region. 

The total oxygen mass is independent of the X-ray lu- 
minosity, while it increases by a factor ~2-3 with increasing 
Lfuv due to enhanced photodissociation of H2O and OH in 
the outer disk. Since the line is optically thick, such changes 
in mass ~ 2-3 do not affect the total line flux. 

The oxygen emitting region is generally pushed fur- 
ther out with increasing Lx, as the disk is overall warmer 
(Paperl). This is also reflected in the FWHM of the 
[0 1] line, being narrower for high Lx- However, high spec- 
tral resolution, such as SOFIA/GREAT is needed (R ~ 
75000) to observe this effect. 
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Fig. 9: Our [0 1] 63 /mi results compared with other models. 



In Fig. [HI we show the results from our model grid for 
[0 1] 63 /im together with the fluxes found b y other authors. 
Our r esults are in good agreement with iMeiierink et al.l 
(|2008T ): their flux at L x = 2 x 10 30 erg s -1 matches our values 
for low Lfuv models. Their prediction f or Lx = 2x 10 29 erg 
s _1 is lower than ours, but as noted in lAresu et all (I201lh 
this is due to the absence of FUV in their X-ray irradiated 
model. They predict a slightly higher flux for their highest 
X-ray luminosity (Lx = 2xl0 31 erg s^ 1 ). 

The lErcolano et "all (120081 ) value (log F ~ -16.41 
W/m 2 for Lx— 2xl0 30 erg s _1 and no FUV) is in very 
good agreement with our predictions for low Lfuv- 
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Fig. 8: Correlation of water lines with the [O i] 63 /im lines. From left to right: 23 /im, 63 /xmand 179 /im. Color coded 
is Lx, different symbols indicates different FUV luminosities. The linear fit for the UV only model and the UV+X-rays 
model are drawn with a black and grey line respectively. 



Our values are somewhat h igher compared with the ones 
of iGorti fc Ho llcnbach (2008|). In their models, they con- 
sider L FU v= 5 x 10 31 ergs _1 and L x ^ 2 x 10 30 ergs _1 (fidu- 
cial model A) . They also calculate a model without X-rays 
(model C) and a model D with 10 times higher Lfuv- They 
consider different absolute values for LpuVj but we repro- 
duce the same [O i] flux variation for models where Lfuv is 
scaled by a factor 10. Increasing the FUV luminosity by an 
order of magnitude (from model A to model D) causes the 
flux to increase by a factor ~ 5, which is the same factor we 
find for such an FUV luminosity variation. The discrepancy 
in the absolute fluxes might arise from further differences 
in the model, e.g. their grain size distribution goes from 50 
A to 20 /Ltm. 

The dust parameters we explored in the grid cause a 
spread in the [Oijflux of less than a factor 3. The same 
holds for the range of surface density power law exponent 
explored here (i.e. 1, 1.5). Their impact on the [Oi]63 
/im line flux is marginal for our grid, but should defi- 
nitely be taken into account if the goal is to perform multi- 
wavelength fitting of an individual source. In that case 
Lx and Z/fuv are no free parameters, and the secondary 
disk parameters play the dominant role. 

4.2. [C II] at 157 /im 

The FUV luminosity controls the line emission, mainly via 
carbon ionization. Thermal effects impact the line similarly 
to what happens for [Ol], but to a lesser extent. The inner 
and outer radius of the C + emitting region respond to the 
FUV radiation. 

Our fluxes are hi gher than the value s calc ulated by 
iMeiierink et al.l (|2008f) and lErcolano et all (|2008l ). because 
these authors do not consider FUV radiation in their mod- 
els. X-rays alone are not able to sustain enough C + produc- 
tion for the flux to be on the observable limit of Herschel. 
However, in many of our models presented here, we overpre- 
dict the C + line: only few sources in Taurus were observed 
with fluxes above 10~ 17 W m~ 2 (Howard et al., 2012, in 
prep.). One reason for this discrepancy might be the shape 
of the FU V spectrum. As no ted already by iBergin et al.l 
(|2003h and iFoeel et~aT1 (|201l[) . the spectral shape in this 
range is far from being the simple power law that we 
adopted for this study. Many emission lines from highly 
ionized species, especially a very strong Lya line (10.1 eV, 



1215 A), contribute to the total £fuv- The lines in a typi- 
cal T Tauri FUV spectrum can carry up to 75% of the total 
flux for energies lower than the carbon ionization potential 
(11.2 eV, ~ 1090 A). The C+ mass for a given L FU v, de- 
pends on how much flux in the FUV range is carried by 
these lines longward of 1090 A. 

In our case of power law continuum, we likely over- 
estimate, for a given LfuVj the continuum flux between 
11.2 and 13.6 eV, and this might lead to an overproduction 
of C + . X-rays do not contribute significantly to the C + 
production, because upon X-ray ionization, C + likely loses 
more than one electron through the Auger effect. X-rays 
also destroy C + via direct ionization, which has roughly the 
same rate as neutral carbon X-ray ionization. Subsequent 
recombination channels form C + , but also neutral carbon. 
X-rays do not sustain the same level of carbon ionization 
as FUV photons (see also the column densities in Paperl). 

4.3. Neon 

[Neil] at 12.8 /im probes the innermost high temperature 
conditions in the top layers of our protoplanetary disk mod- 
els. The X-ray photons are necessary to produce ionized 
neon and heat the gas there. The temperatures needed to 
excite the transition via electron collisions (T cx ~ 1000 K) 
can only be achieved in the innermost ~ 50 AU in a X-ray 
(or EUV) irradiated disk. 

Our models predict a correlation between [Ne n] and X- 
rays; all other parameters explored in this paper contribute 
marginally to the formation of the line. The slopes in the 
fits for M Nc + and -F[Neii] vs -^x are very similar confirming 
that the line is optically thin. The higher the X-ray flux, 
the more extended is the emitting region. This causes the 
FWHM of the line to anticorrelate with X-rays, changing 
from ~ 25 to ~ 15 km s^ 1 when Lx goes from 10 29 to 10 31 
erg s _1 . 

The line i s mainly observed in sources that drive out- 
flows or jets (Sa cco et al.l . l2012t iBaldovin-Saavedra et all 
I2012T ). theoretical calculations have shown that the line 
can be indeed excited in such outflows, providing line 
fluxes and profiles in a ccordance with the observation 
(lErcolano fe Owed . [2010h . 

A line probably em itted purely from a disk was observed 
bv lNaiita etaLI (120091) toward GM Aur (L x = 1.2 xlO 30 erg 
s~ 1 . lGiidel et al.l (|2010t 0. The authors note the line is cen- 
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Fig. 10: Line fluxes for all the tracers discussed in this paper, for different inclinations of the disk: i — (face-on), 45, 60, 
75, 90. 



tred at the stellar velocity and has a FWHM of 14 km/s. 
Our model predictions are consistent with this value (Fig. 
[5]), and confirm that such FWHM values can be readily 
reached. The same line shows a redshifted contribution at 
~ 50 km/s. Interestingly this could be either a contribution 
from a jet component (t hough no known j ets/outflow are 
associated with GM Aur. lGudel et al.l (|2010h and references 
therein) or from the inner rim of the disk (the authors sug- 
gest 0.1 to 0.5 AU). The blushifted counter part was not 
observed because of atmosphere issues and its presence or 
absence can thus not be assessed. In Fig. |S1 we show our 
prediction for a generic model, where the disk inner radius 
is 0.5 AU. Though we are not able to reproduce the high 
velocity shift of the second component, we clearly predict 
the presence of two components in the line profile. The 
high velocity component originates close to the inner rim 
(25 km/s correspond to ~ 1-2 AU in our model), while the 
second lower velocity component originates between 2 and 
5 AU. 



In this case, the emission is in LTE since the electron 
density in the inner rim exceeds the critical density for this 
line. If the high velocity component is confirmed to be sym- 
metric, it could present an interesting diagnostic tool for 
studying the geometry (height, position, shape) of the in- 
ner rim. 



In our models the correlation between [Ne n] and Lx is 
clear, although we rely on one set of fiducial disk param- 
eters and the X-ray spectrum is scaled to match the re- 
quired Lx, leaving its shape unvaried. Observations include 
sources that have different disk properties, which can af- 
fect [N e II] as well. This was investigated by ISchisano et al.l 
([2010D . They could reproduce the 1 dex scattering in LfNeii] 
observed by iGtidel et alJ (|2010h by varying the hardness 
of their X-ray spectrum and the flaring angle of the disk 
model. 
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4.4. Ar+ 

Our [Aril] 6.9 /im results are very similar to the 
[Ne n] results described above. The [Ar n] line traces exactly 
the same region, since the chemistry and excitation of these 
two species are almost identical. We need to stress here, 
that the collisional rates used to calculate the [Ar n] line 
fluxes are a factor ^10 higher than the collisional rates for 
Nc + . Given the same basic chemistry for these two species, 
this is the main reason for the discrepancy in the calculated 
fluxes. Nevertheless, we are interested in the qualitative be- 
haviour of this tracer, which appears to be an interesting 

complem ent to [Nell], 

Also iHollenbach & Gortil (|2009t ) give prediction for 
this tracer concluding that the total fluxes of [Ne n] and 
[Ar n] are comparable. A study of Spitzer IRS archival spec- 
tra was able to detect this line in protoplanetary disks 
(Szulagyi et al., 2012, subm.). More observation are needed 
to confirm the diagnostic of this tracer. Our estimates for 
the ionized ar gon flux are higher than th ose predicted by 
other groups l|Gorti fc Hollenbachl . 12008ft . This might be 
due to the different recipe used for the collisional rates. 

4.5. Influence of inclination 

So far we restricted the discussion to i = 45°, however to 
estimate the flux variation of the tracers discussed in this 
paper, we selected a sub-set of 20 models, in which we vary 
the X-ray and FUV luminosities. The other parameters are 
fixed: a min = 0.1, a pow = 1.5 and e= 1.5. For these 20 mod- 
els, we run the line radiative transfer for four more different 
inclinations (0° face-on, 60°, 75° and 90°). In Fig. [TO] we 
plot the line fluxes versus Lx, color coded are the different 
inclinations at which the disk is seen by the observer. 

The [O i] flux (Fig. [TU1 bottom left panel) is marginally 
affected (0.1-0.2 dex) for Lx < 10 30 erg s _1 and inclinations 
lower than 75°. The flux naturally drops (factor 5-10) for 
an edge-on disk. When Lx > 10 30 the scatter introduced 
by models with i different than 45 is ~ 1 dex. Even con- 
sidering these extreme inclinations, there are no models 
with [Oi] fluxes below 10~ 17 W/m 2 ; hence our threshold 
for [Oi]is very robust. 

The ionized carbon flux and water flux at 179 /imare 
only marginally affected by different disk inclinations, as it 
is expected for outer disk tracers. The [Ne n] flux (as well as 
[Ar n] , not shown) do not depend on the inclination when 
i < 75°, while the flux drops by ~ 1 dex for i = 75° and 
by more than 2 dex for an edge-on disk. Sources with incli- 
nation higher than 60° are then likely to contribute to the 
scatter in the [ Nell] line observed by iGiidel et all (|2010t ) 
and modelled bv lSchisano et al.l (|2010t ). 

The water line emission at 23 /im and 63 /im is sensitive 
to the disk inclination, because the lines are emitted in 
the inner disk. The scatter is of the order of 0.3 dex or 
the face-on and i = 60° models, and ~ 1 dex for i = 90°. 
The inclination is then an important parameter to take into 
account when the correlation between these two water lines 
and [O i] is evaluated from a random sample of sources. 

4.6. Future 

The results of this grid provide the tool to understand 
observations of a larger sample of disks. In the context 
of the GASPS open time key project with Herschel (Gas 



in Protoplanetary disks, P.I. Bill Dent), observations of 
[C n] and [O i] are available for more than 100 sources in 
Taurus (Dent et al., in preparation). The [O i] fluxes ob- 
tained in the GASPS sample are in qualitative and quan- 
titative good agreement with the predictions made in this 
work (Howard et al., in preparation). A detailed compari- 
son with the predictions made in this work is the topic of 
a future paper. 

The Spitzer satel l ite ob served the [Ne n] line in many 
sources (jGiidel et al.l . I2010D . Due to the low spectral reso- 
lution, it is difficult to assess the location of the emitting 
region. Recent observations with high spectral resolution 
instruments provide a chance to disentangle the line origin. 
Currently there are too few sources to carry out a statisti- 
cal study of the disk emission versus jet emission. Clearly 
more high spectral resolution ground based observations 
are needed. 

A range of water lines with different excitation en- 
ergy probe different disk regions, ranging from hot water 
in the inner rim to the outermost cold water that orig- 
inates from photodesorption and ion-molecule chemistry. 
The correlation of water lines with the [OI] fine struc- 
ture line at 63 /xm could give information on the disk 
structu re and the excitation condition s in the disk atmo- 
sphere. iRiviere^Marichalar^r^l] (|2012l ). for example, shows 
how the correlation between the water and the oxygen 
line at 63 /imhas a slope lower than 1 as predicted here. 
Dedicated modelling for t his and other water observations 
(|Pontoppidan et al.Ll2010f ) will also be addressed in a future 
paper. 

5. Conclusions 

As shown in paper I, changes in the disk structure due 
to different luminosity ratios Lx/Lpuv affect the optical 
thickness and height of the inner rim as well as the posi- 
tion of what we define as the "second bump" . Hence, models 
differing only in the X-ray luminosity have different optical 
depth structure. This modifies the shadowing effects on the 
outer region of the disk, and impacts directly the thermo- 
chemical conditions there, hence changing the physical con- 
ditions that drive the line emission of different tracers. 

Here, we presented an extensive description of fine- 
structure emission lines of oxygen, ionized carbon, ionized 
neon, argon and water. 

— The [O i] 63 (im line emission is optically thick, and 
probes thermal conditions of the gas above the molecu- 
lar layers. The line flux increases with temperature. For 
Lx < 10 30 erg s _1 , FUV dominates the thermal bal- 
ance, above that, Lx> 10 30 erg s _1 , X-rays dominate. 

— [Cn]157 /imis mainly driven by Lpuvvia carbon 
ionization. X-rays affect the line flux to a lesser extent 
through Coulomb heating in the [C il] emitting region. 
The detailed emission line spectrum and continuum 
flux levels of the FUV spectrum are important in order 
to explain observed [C n] line fluxes. 

— The [Nell] emission from our static modeled disk 
atmosphere correlates with X-rays. The line probes 
thermal conditions in the upper layers of the disk 
in the region between 1 and 10 AU. Our line profile 
predictions can be tested with high spectral resolution 
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ground based observation. Especially the presence 
of a high velocity double-peaked component in the 
line profile, next to a low velocity double-peaked 
component, could give informations about the position 
of the gas inner rim independent of dust observations. 
This line profile is meant to be compared to sources 
for which disk emission has been established, as 
it cannot reproduce line profiles arising from gas in 
radial motion, such as in photoevaporative flows or jets. 

— [Aril] traces the same disk regions as [Neil]. Recent 
work (Szulagy et al., 2012, subm.) confirm the detection 
of this line in protoplanetary disks in few objects. 

— Water emits over a wide range of radii, from the inner to 
the outer disk. Water line fluxes behave differently de- 
pending on the disk region they arise from. The diverse 
correlations with [O i] 63 fxm can confirm the validity of 
the overall chemical and thermal structure of disk atmo- 
spheres as modelled with thermo-chemical codes such as 
ProDiMo. 
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Fig. A.l: Water emitting regions for a subsample of models with a m j n = 0.1 /im, a pow — 3.5 and e = 1.5. Yellow, cyan 
and green mark respectively the emitting region of the 0-H2O rotational transitions at 23.8, 63.3 and 179 /im. Blue is 
the emitting region of [0 1] 63 fini. In grey scale we plot the water density. From left to right are models with increasing 
Lpuvfrom 10 29 to 10 32 erg s _1 , from upper to lower panel Lx increases from to 10 32 erg s _1 . 
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